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At high temperatures SiC and SLN4 react with water vapor to form a silica 
scale. Silica scales also react with water vapor to form a volatile Si(OH)4 
species. These simultaneous reactions, one forming silica and the other 
removing silica, are described by paralinear kinetics. A steady state, in 
which these reactions occur at the same rate, is eventually achieved. 
After steady state is achieved, the oxide found on the surface is a constant 
thickness and recession of the underlying material occurs at a linear rate. 
The steady state oxide thickness, the time to achieve steady state, and the 
steady state recession rate can all be described in terms of the rate 
constants for the oxidation and volatilization reactions. In addition, the 
oxide thickness, the time to achieve steady state, and the recession rate can 
also be determined from parameters that describe a water vapor-containing 
environment. Accordingly, maps have been developed to show these 
steady state conditions as a function of reaction rate constants, pressure, 
and gas velocity. These maps can be used to predict the behavior of silica- 
formers in water-vapor containing environments such as combustion 
environments. Finally, these maps are used to explore the limits of the 
paralinear oxidation model for SiC and SLN4. 


INTRODUCTION 

SiC and SLN4 are proposed for long-term applications in combustion environments. 
Some of these applications include combustor liners and turbine vanes for both 
propulsion and power generation. The combustion environment for hydrocarbon/air 
systems contains about 10% water vapor, independent of hydrocarbon type and fuel-to- 
air ratio (1). The reactions of SiC and SLN4 with water vapor are therefore a concern. 
While much of the information in this paper is generic for any environment containing 
water vapor, this paper will emphasize the combustion environment in a gas turbine since 
the work summarized here was driven by development for these applications. 

In water vapor -containing environments SiC and SLN4 undergo both an oxidation and a 
volatilization reaction, as shown below for SiC: 

SiC + 3H 2 0(g) = Si0 2 + 3H 2 (g) + CO(g) [ 1 ] 

Si0 2 + 2H 2 0(g) = Si(OH) 4 (g) [2] 
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The kinetics of the oxidation reaction are described by the parabolic rate constant for 
oxide formation, k p , whereas the kinetics of the volatilization reaction are described by 
the linear rate constant for oxide volatilization, k|. These reactions occur simultaneously 
and are described by paralinear kinetics. The paralinear kinetic model has been 
developed for simultaneous oxidation and volatilization of chromia-formers by Tedmon 
(2) and is directly applicable to the oxidation of silica-formers in water vapor: 


dx 

d7 



[3] 


where x is oxide thickness and t is time. At long times or high volatility rates a steady 
state is achieved in which oxide is formed at the same rate it is volatilized. SiC and Si 3 N 4 
undergo linear recession, given by the rate yL, when this steady state is achieved which 
is directly related to the volatility rate of the oxide, kj. Dimensional changes due to 
paralinear oxidation are shown as a function of time for SiC >2 growth on SiC in Figure 1. 
At steady state, a limiting oxide thickness, xl, is achieved which is given by (2): 


The subscript L here and in all future uses refers to the steady state limit. The paralinear 
oxidation kinetics of SiC and Si}N 4 in water vapor and combustion environments have 
been studied under a wide variety of conditions (3-10), and the model has been shown to 
be valid in a majority of cases. It is therefore timely to extend the description of steady 
state oxidation/volatilization kinetics and compare the theory to observations obtained for 
SiC and SIjN 4 over a wide range of experimental conditions. This work was begun in a 
previous paper (11) but is more completely addressed herein. 

Specifically, the aims of this paper are as follows. First, the concept of the time to reach 
steady is state is developed to extend Tedmon’ s treatment of paralinear kinetics. Second, 
general relationships for limiting oxide thickness, the time to reach steady state and 
steady state recession rate are developed in terms of the total pressure, gas velocity and 
temperature in the combustion environment. Third, a series of constant temperature maps 
are presented which graphically show the variation of the steady state oxide thickness, 
time to reach steady state, and steady state recession rate with both rate constants and 
combustion environment parameters. Experimental results are then compared to 
theoretical predictions. Finally, these maps are used as a framework for discussing the 
limitations of the paralinear oxidation model for SiC and ShlNU in combustion 
environments. 


DETERMINATION OF THE TIME TO ACHIEVE STEADY STATE 

In Tedmon’s treatment of paralinear kinetics, the steady state oxide thickness, as well as 
the steady state recession rate, are described mathematically (2). While the time to reach 
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95% of the limiting oxide thickness is shown graphically for a real chromia former as a 
function of temperature, a mathematical expression was not developed. In describing 
paralinear kinetics for a system, the time to reach steady state is essential for predicting 
when the steady state oxide thickness is formed, as well as when a simple linear rate 
describes the recession of the underlying material. In order to determine the change in 
oxide thickness with time, the integrated form of Equation 3 is needed. This expression 
is: 


t = ^ 


2kf 


2k, 


x-ln 
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yj 


[5] 


Upon inspection of this equation as well as Figure 1, it can be seen that the oxide 
thickness, x, approaches the limiting thickness, xl, given by Equation 4, asymptotically 
with time. An expression for the time at which the limiting oxide thickness is achieved 
cannot be defined. Instead, as is typical for this kind of a function, a time constant can be 
defined which is that time at which the oxide thickness given by parabolic kinetics (valid 
at short times) is equal to the oxide thickness found at steady state 
oxidation/volatilization (valid at long times). This time constant then defines the 
transition time at which the simple parabolic oxidation model becomes less accurate than 
the simple steady state oxidation/volatilization model. The graphical determination of the 
time constant for establishing a limiting oxide thickness, t Cl0 , is shown in Figure 2 and is 
given by: 


* 0,0 — ‘ 


4(k,r 


[ 6 ] 


In Figure 2, both the oxide thickness and time have been normalized by the values Xl and 
t Cj0 so this figure is general for any rate constants, k p and k|. Independent of rate 
constants, the time constant represents that time at which the oxide thickness is 70% of 
the limiting oxide thickness, xl. In order to define a time to form an oxide thickness that 
is closer to the limiting oxide thickness, some multiple of the time constant t co must be 
used. For example, for the purposes of this paper it was decided that a multiple of 2 x t c>0 , 
which represents the time to reach 84% of the limiting oxide thickness, is a better 
description of the time to reach nearly steady state oxidation/volatilization. This time has 
been defined as the limiting time, tL, i.e.: 


t L = 



[7] 


A table is included in Figure 2 that shows the multiples of the time constant that 
correspond to the time to form a certain percentage of the limiting oxide thickness. 
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Similarly, a time constant can be formulated to describe the time to achieve steady state 
recession. In this case, the time constant is defined as that time at which the expression 
for parabolic recession (valid at short times) intersects with the expression for linear 
recession (valid at long times). This time constant then defines the transition time at 
which the simple parabolic recession model becomes less accurate than the simple steady 
state recession model. The graphical determination of the time constant for establishing a 
steady state recession rate, t c , r , is shown in Figure 3 and is given by: 


t 


c,r - 



[B] 


Here k’ p and k’i represent the rate constants for parabolic and linear recession of SiC. 
(These are different from the rate constants for parabolic oxidation and linear 
volatilization of SiC >2 used in Equations 3 through 7. The units k’ p and k’i are in |i 2 /h and 
|i/h of SiC or SLN 4 consumed rather than p 2 /h and |i/h of SiCL formed or volatilized as in 
Equations 3 through 7. The rate constants are related by the molecular weights and 
densities of silica and the underlying material as shown in the next section.) Here, both 
the recession depth and time have been normalized by y n and t c , r . The quantity y n is given 
by k p ’/ki’ and is somewhat arbitrarily chosen. Figure 3 is then general for any rate 
constants k p ’ and k|’. Figure 3 contains a table that shows the multiples of the time 
constant that correspond to the time at which the recession rate corresponds to a 
percentage of the steady state linear recession rate. At one time constant, for example, 
the observed recession rate for paralinear oxidation/volatilization is 85% of that which 
would be found at steady state oxidation/volatilization. For the purposes of this paper 
one time constant has been chosen as an adequate description of the time to reach nearly 
steady state linear recession rates. 


PARAMETRIC EXPRESSIONS FOR OXIDATION AND VOLATILIZATION 

In the previous sections the steady state oxide thickness, the time constant to reach steady 
state, and the recession rate have been shown to depend on the oxidation rate as well as 
the volatilization rate. It is also useful to understand how the oxidation/volatilization 
processes in turn depend on more fundamental parameters (pressure, water vapor partial 
pressure, gas velocity and temperature) that describe any environment, rather than on 
empirical kinetic rate constants that may or may not be known. In order to understand 
how the limiting oxide thickness, time constant and recession rate depend on these 
parameters, the parametric expressions for the rate constants must be known. These are 
discussed in the following sections. 

Parametric Expression for the Parabolic Oxidation Rate Constant 

The parabolic rate constant describes the kinetics of Equation 1 for moderately thick 
scales (diffusion limited growth). This rate constant varies with the water vapor partial 
pressure and temperature but is independent of gas velocity. The temperature 
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dependence will be discussed later. The water vapor dependence can be described using 
a power law: 


where n is the power law exponent. This power law should be identical for all silica - 
formers since parabolic oxide growth is rate limited by transport of oxidant through the 
silica scale. Deal and Grove (12) observed a power law exponent of unity for the 
parabolic oxidation of silicon by water vapor, indicating that the oxidation reaction is 
limited by the diffusion of molecular water, rather than a charged species, through the 
silica scale. Studies on the oxidation of SiC in water vapor found power law exponents 
varying between 0.76 and 0.91 (13,14). These parameters are summarized in Table 1. 
Studies on the diffusivity of water vapor in silica also find that molecular water vapor is 
the diffusing species (15). For the purposes of this study a power law exponent of one 
will be used to describe the water vapor partial pressure dependence of the parabolic 
oxidation rate constant for SiC and Si.^lNU. Using results obtained for SiC oxidation at 
1316 °C (13) the following expression is used in this paper to extrapolate oxidation rates 
in water vapor for silica formers: 


This temperature was chosen as a reference temperature since this is the only temperature 
where a complete data set is available for both rate constants. 

Parametric Expression for the Linear Volatilization Rate Constant 

The linear volatilization rate constant describes the kinetics of Equation 2. This rate 
constant varies with the water vapor partial pressure, total pressure, temperature and gas 
velocity. The temperature dependence will be discussed later. Under combustion 
conditions, the volatilization of silica to form Si(OH )4 is limited by transport of Si(OH )4 
through a laminar gaseous boundary layer. This rate is given by the following expression 
for a flat plate geometry: 


where Re is the dimensionless Reynold’s number, Sc is the dimensionless Schmidt 
number, D is the interdiffusion coefficient of Si(OH )4 in the boundary layer gas, p v is the 
equilibrium concentration of the volatile Si(OH) 4 , and L is a characteristic length. This 
equation can be further expanded: 


k P’1316°C -0 - 44P h 2 O ’ [P' Si °2 /h ] 


[ 10 ] 



[ 11 ] 



[ 12 ] 
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where v is the gas velocity, r\ is the gas viscosity, and p is the concentration of the 
boundary layer gas. This expression can then be simplified to the combustion 
environment parameters of pressure and velocity: 


kl 


-.1/2 p 

v ^Si(OH) 4 
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total 


[13] 
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and knowing that Psi(OH) 4 i s proportional to Q from Equation 2, and Ph-,o * s e q ua l 
to about 0. 1 P t otai for combustion environments as discussed in the introduction, then: 


kl 
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1/2 


p3/2 

total 


[14] 


Remarkable agreement between the empirical results and boundary layer transport model 
has been found (7,8). The measured dependencies of volatilization are also summarized 
in Table 1. Using results obtained for SiC oxidation/volatilization at 1316 °C (7) the 
following expression is used in this paper to extrapolate oxidation rates in water vapor for 
silica formers: 


k l’1316°C 


0.36 



p l/2 

total 


. 1/2 


, [|i Si0 2 /h] 


[15] 


The pressure dependence will vary, of course, in conditions where other volatile species 
such as SiO(g), SiO(OH) 2 (g), or SiO(OH)(g) are more stable (8). 


Parametric dependence of limiting oxide thickness, time constant, and recession rate 

It is now straightforward to substitute the parametric dependence of the oxidation and 
volatilization rate constants (Equations 10 and 15) into the previous expressions for xl 
and tL (Equations 4 and 7) to derive the following expressions: 

pl/2 

*L.|3I6"C =0.615 ,[[lSi0 2 ] [16] 

P H ,0 v 1 ' 2 

And for hydrocarbon combustion environments where Ph-,o > s 0-1 Ptotai. the following 
general relationship holds true: 


*L“ (Petal Vi " 1 ' 2 


[17] 
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The steady state oxide thickness therefore decreases with both pressure and gas velocity. 
Similarly, 


t L>1316°C “ total - , [h] [18] 

v 

H z O v 

And again, for hydrocarbon combustion environments where P^q > s 0- 1 Ptotai, the 
following general relationship holds true: 


tL 



[19] 


The time to achieve steady state also decreases with both pressure and gas velocity. 


As stated earlier, the steady state recession rate, y L , is directly proportional to the 
volatility rate of silica, k|, and is given by: 


. _ k p(Si0 2 ) MW(SiC) 
YL_ ' MW(Si0 2 ) p(SiC) 


[ 20 ] 


Here p is the density and MW is the molecular weight of the indicated species. The exact 
semi-empirical relationship for yL is: 


P 2 v 1/2 

y L .1316°C = 0-18 J1 ^ ,[pSiC/h] [21] 

^total 

And again, for hydrocarbon combustion environments where P^o ' s 0-1 Ptotai, the 
following general relationship holds true: 

y L ocP 3/2 v 1/2 [22] 

For easy reference the general expressions for the steady state conditions are summarized 
in Table 2. 


CASE STUDIES FOR PARALINEAR OXIDATION OF SiC AND Si 3 N 4 

Up to this point the theoretical expressions for paralinear oxidation and volatilization; the 
parametric dependence of the steady state oxide thickness, time constant, and recession 
rate; as well as semi-empirical expressions at 1316 °C have been developed. It is now 
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appropriate to introduce a number of studies that have been conducted over a wide range 
of conditions that will be compared to the theory in the remainder of this paper. These 
case studies include three furnace studies (3-5) as well as three studies in actual 
combustion environments (6,7,10). The test parameters for each of these studies are 
listed in Table 3. The results obtained in a high-pressure burner rig at NASA Glenn (7) 
were used to establish the empirical constants for Equations 16, 18, and 21, although the 
parametric dependence in these equations is derived from theory (and verified by the 
experimental observations). This is not a complete list of studies that describe the 
paralinear oxidation/volatilization of silica-forming materials. Other excellent studies 
can be found in References 17 and 18. The conditions for these studies are similar to 
those in Reference 7, thus they were not included in the case studies. 


LIMITING OXIDE THICKNESS AND TIME CONSTANT MAPS IN TERMS OF 

RATE CONSTANTS 

The concept of an oxide thickness map can first be developed in terms of, k p , the 
parabolic oxidation rate constant and, k|, the linear volatilization rate constant. This is 
shown in Figure 4, calculated for SiC at 1316°C. Each diagonal line is a line of constant 
limiting oxide thickness. Thus for a given oxidation rate and a given volatilization rate, 
the steady state oxide thickness at 1316°C can be picked directly from this map. A 
number of points are shown on this map that plot limiting oxide thickness for each case 
study. Steady state oxide thicknesses were calculated using Equation 4. Experimentally 
determined rate constants were not available in all cases, so to be consistent, calculated 
rate constants were used in all cases. Rate constants were calculated using Equations 10 
and 15 and the experimental parameters in Table 3. Calculated rate constants are 
compared to those experimentally measured in Table 4. In most cases there is very good 
agreement between calculated and observed rate constants. This comparison will be 
discussed in greater detail in a following section. 

Similarly, a steady state time constant map as a function of rate constants is shown in 
Figure 5. Here again, rate constants for each case study are calculated from Equations 10 
and 15 based on parameters in Table 3 and are used in Equation 7 to plot steady state 
time constants. Each diagonal line is a line of constant time to reach steady state. Thus 
for a given oxidation rate and a given volatilization rate, the steady state time constant for 
SiC at 1 3 1 6°C can be picked directly from this map. 

It is unnecessary to construct a map of steady state recession as a function of rate 
constants since it depends on the linear volatilization rate of silica alone. 

These two maps, Figures 4 and 5, demonstrate a number of concepts. First, these maps 
allow easy visualization of the competition between oxide growth and oxide 
volatilization and the effects on the expected surface oxide thickness as well as the rate at 
which the steady state paralinear process is approached. Second, the case studies really 
cover a wide range of rate constants and this will allow limits for the paralinear model to 
be explored in detail. Third, the time constant is much more sensitive to changes in rate 
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constant than the limiting oxide thickness. Time constants vary from 2 minutes to 4xl0 3 
hours. This is seven orders of magnitude variation in time to reach steady state. The 
limiting oxide thickness varies by three orders of magnitude in this same rate constant 
space. Finally, these maps show that it may be more useful to think of how the limiting 
oxide thickness, time constant, and steady state recession rate vary with parameters that 
describe the combustion environment. Experimentally determined rate constants are not 
always available. In addition, these maps do not help to visualize how the paralinear 
process is affected by changes in pressure or gas velocity. The development of maps as a 
function of pressure and velocity is the topic of the next section. 


LIMITING OXIDE THICKNESS, TIME CONSTANT, AND RECESSION RATE 
MAPS IN TERMS OF TOTAL PRESSURE AND GAS VELOCITY 

The concept of oxide thickness maps, first developed in terms of the oxidation and 
volatilization rate constants, can now be extended to the parametric expression developed 
in Equation 16 above. Figure 6 shows the limiting oxide thickness contours calculated for 
SiC at 1316°C as a function of total pressure and gas velocity. Each diagonal line is a 
line of constant limiting oxide thickness. Thus for a given total pressure and a given gas 
velocity, the steady state oxide thickness for SiC at 1316°C can be picked directly from 
this map. The points shown on this map plot actual gas pressures and velocities for five 
case studies for SiC or Si.-*N 4 exposures in experimental or combustion environments with 
the steady state oxide thickness calculated using Equation 16. (The case study for 
laboratory furnace 2 is omitted on these maps since the simplification that Ph 2 o = 0-1 

Ptotai can not be made for these test conditions. In addition, it was assumed the water 
vapor partial pressure in the high pressure furnace was 1 atm so that the assumption that 
P H , 0 = 0. 1 Ptotai could be made for this case.) 

Similarly, maps for time to reach steady state and steady state recession rate as a function 
of total pressure and gas velocity calculated for SiC at 1316°C are shown in Figures 7 
and 8, respectively. These maps were calculated using Equations 18 and 21. The points 
on these maps again represent the actual total pressure and gas velocity for the five case 
studies. 

While these maps (Figures 6 through 8) give less information about the competing silica 
formation and removal processes in the paralinear mechanism than Figures 4 and 5, they 
do allow easy visualization of the effects of changes in total pressure and gas velocity on 
the expected surface oxide thickness, the rate at which the steady state paralinear process 
is approached, and the steady state recession rate of structural material. Again, these 
maps demonstrate that the steady state time constant is most sensitive (closely spaced 
contours) to changes in pressure and gas velocity, the recession rate has intermediate 
sensitivity, and the limiting oxide thickness is least sensitive to changes in combustion 
gas parameters. 
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COMPARISON OF EXPERIMENTAL RESULTS TO THE PARALINEAR MODEL 
Rate constants 


As mentioned earlier, the measured and calculated rate constants are compared in Table 
4. In general there is very good agreement. The rate constants measured for SiC in all 
cases are within a factor of two of the calculations based on the semi-empirical 
expression from the high pressure burner rig results at 1316°C, Equations 10 and 15. 
Several points should be made. First the experimental rates were obtained at various 
temperatures as reported in Table 3, whereas all calculated rates are for 1316°C. 
Nevertheless the agreement is quite good indicating the empirically derived constants in 
Equations 10 and 15, as well as the parametric dependence, work well over a wide range 
of conditions. Second, the mechanism for oxidation in the high pressure furnace can not 
be described using a parabolic rate constant due to the porous oxide morphology formed 
(see Figure 9c), so this case is not well represented by paralinear kinetics. This case will 
be discussed in more detail in the following sections. The good agreement between 
calculated and measured rate constants in the majority of the cases gives confidence that 
the derivations of limiting oxide thickness, time constant, and steady state recession rate 
that are based on these rate constants can be used to make accurate predictions of 
paralinear oxidation/volatilization behavior. 

Limiting Oxide Thickness 

Figure 9 shows cross-sections of oxides formed on SiC or Si.^ for each case study 
shown in Table 3. Table 5 summarizes these findings as well as the predicted limiting 
oxide thickness for each case study. 

The oxide grown on SiC in laboratory furnace 1 (Figure 9a) is about 2 |im thick. This is 
much less than the calculated limiting oxide thickness of 91 pm, however, the exposure 
time is only 100 hours which is only a miniscule fraction of the time required to achieve 
steady state. The oxidation kinetics in this case were well modeled by parabolic kinetics 
alone since the volatility of silica is negligible. 

The oxide grown on SiC in laboratory furnace 2 (Figure 9b) is, on average, about a factor 
of two less than the predicted limiting oxide thickness. The exposure time was one-third 
the time to achieve steady state. These results are in good agreement with expectations. 
In this case, the full paralinear model is needed to model the observations. Both the 
oxidation rate and the volatility rate are needed to model the observed kinetics. 

Figure 9c shows the oxides grown on both SiC and SpN 4 in the high-pressure furnace 
environment. There are several phenomena that must be considered when comparing the 
observed oxide thickness to the predicted limiting oxide thickness. First, the gas velocity 
in this furnace is nearly stagnant (0.05 cm/sec). Silica volatility will be negligible in this 
case. This can be observed quite readily from Figure 4, where the oxidation rate is quite 
high due to the high water vapor partial pressure, but the volatility rate is quite low due to 
the very low gas velocity. In this case, it would be expected that parabolic kinetics alone 
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would be sufficient to model the observed oxide growth. This brings up the second 
phenomenon that must be considered. From Equation 10, it is shown that the parabolic 
oxidation rate has a strong dependence on the water vapor partial pressure. As the 
oxidation rate increases, the rate of production of gaseous oxidation products also 
increases as can be seen from the oxidation reaction (Equation 1). The high pressures of 
these gases generate porosity in the oxide layer as shown in Figure 9c for SiC. The 
amount of pores in the silica scale has been shown to increase with water vapor partial 
pressure (14). Because the scale is no longer dense, the transport path length of water 
vapor through the silica scale is no longer parabolically increasing. Gas phase transport 
can now occur through the interconnected cracks and pores. The oxide growth rate is 
limited by solid state transport of water vapor through the dense portion of the silica scale 
near the interface which remains at a nearly constant thickness with time (5). The 
oxidation kinetics and corresponding recession rates are therefore linear. In the case of 
the SisN 4 exposed in the high pressure furnace, the dense layer is much thicker 
presumably due to the relatively lower amount of gaseous products released per mole of 
silica formed. In this case, parabolic kinetics were observed. The surface layer of oxide 
is still porous and friable. At long times, spallation of this porous surface layer may limit 
the thickness of the silica rather than volatility. It remains to be demonstrated if there is 
some upper limit of silica thickness at which spallation typically occurs. The important 
point to take away from these exposures is that the underlying material and the resulting 
oxide morphology affect the kinetic rate law applicable to oxide growth. A paralinear 
oxidation/volatilization model is not applicable for the conditions in the high pressure 
furnace. 

Figure 9d and 9e show the oxide morphology found for SiC exposed in a mach 0.3 burner 
rig and a high pressure burner rig, respectively. In both cases, the predicted steady state 
oxide thickness is less than that actually observed. This could be explained by impurity- 
enhanced oxidation. It is known that even small amounts of impurities can increase the 
oxidation rate of SiC by an order-of-magnitude (19). Burner rigs and other real 
combustion environments do contain more impurities than laboratory furnaces. If 
impurity effects caused k p to increase by an order-of-magnitude over the rate given by 
Equation 10, then the limiting oxide thickness would be 6|im in both cases. This is in 
good agreement with observations for oxides formed in both burner rigs. 

Figure 9f shows a cross-section of a S(}N 4 vane exposed in an industrial turbine. No 
oxide is visible on the surface. The surface roughness appears to correspond to the shape 
of individual SbN 4 grains with the grain boundary phase removed. The gas velocity in 
this environment is much higher than found in the other case studies. Again several 
phenomena must be considered when comparing the predicted limiting oxide thickness to 
the lack of oxide actually observed. First, for very thin silica layers, parabolic oxidation 
kinetics is not expected. The chemical reaction rate of the water vapor with the 
underlying material could be expected to control the oxidation rate since transport across 
a thin layer would be rapid. In this case, linear oxidation rates would be observed. 
Secondly, because the volatilization rate is so rapid and the time to reach steady state is 
so short, the oxidation rate is only of interest at very short times for an 
oxidation/volatilization mechanism rendering the first point of little interest. Third, if the 
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volatilization is so rapid that no oxide is found on the surface, the kinetics of recession 
may not even be related to the volatility of silica. Several scenarios can be envisioned. 
The recession could be controlled by the rate of chemical reaction of the oxidant with the 
SiC or SfiN 4 to form silica. As soon as oxide is formed it is swept away. Another 
possibility is that a different reaction sequence occurs rather than the reactions in 
Equations 1 and 2. Some kind of “active” volatilization could occur similar to the active 
oxidation mechanism that occurs for silica-formers at low oxygen partial pressures. In 
the case of active oxidation at low oxygen partial pressures, SiC reacts with oxygen to 
form SiO(g) rather than SiC> 2 . An analogous reaction can be envisioned in which SiC 
reacts directly with water vapor to form Si(OH) 4 (g) without the intermediate silica 
formation step. It might be possible to distinguish this kind of reaction by determining 
the temperature dependence of recession. It would be difficult to perform these kinds of 
experiments since these very high gas velocities found in the turbine cannot be readily 
obtained in a laboratory or rig environment. The important point to take away from these 
turbine exposures is that a paralinear oxidation/volatilization model is not applicable for 
very high velocity conditions where the silica volatilization rate exceeds the silica 
formation rate. 

Steady State Recession Rate 

Steady state recession rates calculated from Equation 21 are shown in Table 6. While 
sufficient data are available to compare experimentally observed oxide thicknesses with 
predicted limiting oxide thicknesses, the same is not true for recession rates. Recession 
by paralinear oxidation/volatilization is too small to be measurable in reasonable amounts 
of time for the furnace studies and the mach 0.3 burner rig. Given a one hundred hour 
test and the calculated rates listed in Table 6, recession would only be measurable 
(>25|J.m) for exposures in the high-pressure burner rig and the industrial turbine. 

Material recession has been measured for the materials tested in the high-pressure 
furnace, but as already discussed, this recession can be attributed to consumption of 
material by a linear oxidation mechanism for SiC and parabolic oxidation which may be 
limited by spalling for Si.dvU since volatilization of silica is minimal. 

Equation 21 was determined from the high-pressure burner rig results, so an independent 
comparison cannot be made in this case. However, the good agreement between the 
predicted parametric dependence of the recession rate (Equation 14) and the observed 
recession rate dependence in Equation 2 1 provides some measure of confidence that this 
equation is correct. In addition, the recession measured for silica formers in other 
pressurized burner rig studies (17,18) is in good agreement with Equation 21 for 
moderately different conditions. 

Finally, the calculated rates are in good agreement with the measured rates for the ShN 4 
exposed in the turbine. Because no silica was observed on the surface of the ShN4 vane, 
the paralinear oxidation/volatilization model should not be applicable. However, the 
good agreement of calculated and predicted steady state recession rates may indicate that 
the conditions under which the vane was exposed are just at the limit where the paralinear 
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model may be applied; i.e., at lower volatility rates paralinear kinetics apply and at higher 
volatility rates some of the mechanisms described above in the limiting oxide thickness 
section may be operative. 

While recession rates are difficult to measure when volatilization is only moderately 
rapid, linear weight losses that correspond to the linear recession rates are quite easily 
obtained. The volatilization rate constants shown in Table 4 were obtained by weight 
measurements for cases B, D and E. Measured weight loss and material recession were 
both obtained in the high pressure burner rig and were found to differ by a factor equal to 
the density of SiC, as expected. 


TEMPERATURE DEPENDENCE OF PARALINEAR OXIDATION - 

VOLATILIZATION 

Maps for limiting oxide thickness, time constant, and steady state recession rate as a 
function of rate constants or pressure/velocity space are calculated at a single 
temperature. Different maps must be calculated for each temperature. Up to this point, 
all calculations have been limited to one temperature, 1316°C, and the temperature 
dependence has been neglected. This has been done because oxidation and volatilization 
rate constants measured as a function of pressure and gas velocity for a single material in 
a single rig are not available over a range of temperatures. The available measured 
temperature dependencies of the rate constants are summarized in Table 1. The 
temperature dependence for both rate constants determined in a single study can be found 
only for SiC studied in a laboratory furnace (4) and the uncertainties of the temperature 
dependence in this study are quite large. One point to note is that the enthalpies for 
oxidation and volatilization are similar, so predicting the temperature dependence of the 
limiting oxide thickness and the time to reach steady state, which depend on both k p and 
ki is difficult. However, when good oxidation/volatilization results over a wide range of 
temperatures become available for a given material the following discussion will be 
applicable. 

The temperature dependence of the parabolic rate constant arises from the enthalpy of 
diffusion for molecular water in silica. The values reported for this oxidation enthalpy, 
AH 0X , in the literature are summarized in Table 1 and range from 28 to 156 kJ/mol. The 
temperature dependence of the volatilization reaction arises mainly from the enthalpy of 
Reaction 2. This temperature dependence is found in the p v term in Equations 1 1 and 12 
which arises from the temperature dependence of the Psi(OH) 4 term i n Equation 13. The 

enthalpy for volatilization, AH VO i, varies between 3 and 108 kJ/mol as reported in the 
literature and summarized in Table 1. Thus the temperature dependence of these two 
mechanisms may be similar. 

The temperature dependence of the steady state oxide thickness is given by: 
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X L = 


2k, 


exp 


AH V0 | AH 0X 
RT 


[23] 


where AH is the enthalpy, the subscript vol indicates volatilization, the subscript ox 
indicates oxidation, R is the gas constant, and T is absolute temperature. In reference 4, 
AHvoi = 3±75 kJ/mol and AH 0X = 35±147 kJ/mol. Thus it is highly uncertain which 
enthalpy is larger. If, in fact, AH VO i < AH 0X , then xl increases as T increases. This 
discussion assumes the major volatile species is Si(OH) 4 (g). At high temperatures 
(>1400°C) SiO(OH) 2 (g) or SiO(OH)(g) species may be important (8). These species have 
much higher enthalpies of formation from silica and water vapor. In these cases the 
volatilization reaction would become increasingly important as the temperature was 
increased. The steady state oxide thickness might therefore decrease with temperature. 
There are no oxide thickness measurements available to see how the steady state oxide 
thickness varies with temperature. 

Similarly, the temperature dependence of the time to reach steady state is given by: 




ti =- 


2k; 


exp 


2AH vol - AH 


ox 


RT 


[24] 


where all terms have been previously defined. Again, assuming 2AH VO i < AH 0X based on 
reference 4, then tL will increase as T increases. Clearly, these temperature trends will 
depend on the material of interest and cannot be predicted without better data. 

Finally, the temperature dependence of the steady state recession rate will depend only on 
the enthalpy of volatilization and it is quite straightforward to state that recession rates 
increase exponentially with temperature. 


LIMITS OF THE PARALINEAR MODEL FOR SiC AND Si 3 N 4 OXIDATION - 
VOLATILIZATION IN COMBUSTION ENVIRONMENTS 

The maps developed in Figures 6, 7, and 8 provide a framework for understanding the 
limits of the paralinear model. Several points have been made in the previous discussion 
about the applicability of the paralinear oxidation and volatilization model for silica- 
forming materials in the specific case studies. These points are summarized below and 
in Table 7. A general map, shown in Figure 10, is then developed to show approximate 
pressure and gas velocity limits for applying the paralinear model. Some of these limits 
are based on the long-term application (10 4 to 10 5 hours) of SiC and Si 3 N 4 in a 
combustion environment. Clearly different limits would be chosen for significantly 
different applications. 
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1 . Negligible volatilization/recession 

At low gas velocities and low pressures in a combustion environment, silica volatilization 
and the corresponding steady state recession rate are negligible. Parabolic oxidation 
kinetics are sufficient to model the oxide growth and recession. In Figure 10, a steady 
state recession rate of 1 x 10~ 3 |i/h has been chosen as negligible. For recession rates 
larger than this, the paralinear model should be used. 

2. Long time constant 

At low gas velocities and low pressures in a combustion environment, the time to reach 
steady state oxidation and volatilization can be as long as the expected use life of a 
component. Clearly it is not necessary to apply paralinear kinetics to model the oxidation 
of materials under these conditions. In Figure 10, a time constant of 10 4 hours was 
chosen as a limit for application of the paralinear model. For shorter time constants, the 
paralinear model should be used. 

3. Negligible volatilization, high oxidation rate 

At low gas velocities and high pressures in a combustion environment the volatilization 
rate is low. The oxidation rate is high due to the high water vapor partial pressure. While 
paralinear kinetics can still be used to model the oxidation and volatilization reactions 
under these conditions, the kinetics are dominated by the oxidation reaction and can 
therefore be modeled by the oxide growth kinetics alone. In addition, because oxidation 
rather than volatilization dominates the reaction kinetics, the observed oxidation behavior 
is material dependent. In volatility-dominated conditions (high gas velocity), the 
volatilization of silica appears to be independent of the underlying material or the purity 
of the silica. However, when the kinetics are dominated by high oxidation rates, as they 
are at low gas velocities and high water vapor partial pressures, the oxidation rate and 
mechanism are affected by the underlying material. The purity of the material and the 
product gases generated by the oxidation reaction affect the oxidation mechanism. This 
is evident from the different behavior observed for SiC and ShN 4 in the high pressure 
furnace (Figure 9c). These differences in oxidation mechanism are now discussed. 

3a. Nonprotective scale 

For SiC exposed under low gas velocity, high pressure conditions for 500h, a thick 
porous scale was formed with a 5 micron thick dense silica layer near the oxide/SiC 
interface (5). This dense layer did not grow with time, and the observed kinetics were 
linear. The following mechanism explains these observations. The high oxidation rates 
which occur at high water vapor partial pressures generate large amounts of product 
gases (CO and H 2 ). These gases form bubbles and pores in the silica scale (5, 14). The 
oxidation rate is limited by transport of water vapor through the 5 micron thick dense 
layer of silica near the oxide interface. Thus, a 5 micron steady state oxide thickness has 
been chosen as a limit for the paralinear model in Figure 10. For SiC, only conditions 
which result in a steady state oxide thickness less than 5 microns can be modeled by 
paralinear kinetics (or simple parabolic kinetics). Oxides thicker than 5 microns would 
be full of bubbles and pores and grow by linear kinetics. 
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3b. Thick spalling scale 

For Si}N4 exposed under low gas velocity, high pressure conditions for 500h, 20-25 pm 
thick scales were formed which had very friable surface layers (20). The oxidation 
kinetics for SpN4 oxidized under these conditions were parabolic. The friable surface 
layer indicates that spallation will be a problem for very thick oxides. If spallation occurs 
at some limiting thickness, then kinetics would be parabolic (as a simple case of 
oxidation/volatilization paralinear kinetics) until this oxide thickness is achieved. It is 
possible that as the scale grows thicker, the surface layer spalls off to some average 
thickness. This remains to be verified. If spallation of thick scales occurs, then the 
paralinear oxidation/volatilization model would not be applicable for predicted steady 
state oxide thicknesses greater than the spallation limiting value. In Figure 10, a 
spallation limited oxide thickness of 30 pm has been chosen as a limit for the paralinear 
oxidation/volatilization model. For thicker oxides, a paralinear oxidation/spallation 
model might be a valid. Oxidation/spallation models have been used to describe the 
kinetics of cyclic oxidation of superalloys (21) for many years. A better understanding of 
silica scale growth in high pressure water vapor for very thick scales is needed. 

4. Thin or missing scale 

At high pressures and high gas velocities in a combustion environment silica 
volatilization rates are more rapid than oxidation rates. Very thin scales are formed or 
bare surfaces are observed. The rate limiting mechanism is unknown under these 
conditions and may be any of the following as discussed in more detail above: linear 

oxidation with simultaneous linear volatilization for thin scales, volatilization limited by 
the oxidation rate, or an active volatilization process in which the silicon-based ceramic 
reacts directly to form a volatile hydroxide. In Figure 10, a steady state oxide thickness 
of 0.1 micron has been chosen as the limit for paralinear oxidation/volatilization. For 
steady state oxide thicknesses larger than this, the paralinear model is applicable. 

The combination of all the limits discussed in the four points above maps out a region in 
which the paralinear model for oxidation/volatilization can be applied. In this region, the 
equations summarized in Table 2 are valid. It should be noted that the high-pressure 
boundaries at both very high and very low gas velocities are tentative due to the limited 
number of conditions under which SiC or SpN4 have been tested. The positions of these 
lines are indeterminate, although the slopes should be valid. More testing is needed 
before these boundaries can be definitively established. Many of the propulsion and 
power generation applications for SiC and ShN4, however, do operate in the pressure and 
gas velocity space mapped out by the paralinear regime shown in Figure 10. The mach 
0.3 burner rig and high pressure burner rig both fall within this regime and are therefore 
good tools for screening materials for long-term applications in combustion environments 
where paralinear oxidation/volatilization occurs. 


SUMMARY AND CONCLUSIONS 

Steady state paralinear oxidation and volatilization of SiC and S4N4 in combustion 
environments can be summarized by the relationships given in Table 2. It is important to 
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characterize the oxide morphology in order to understand the steady state recession 
mechanisms and rates. Time constants to achieve steady state are very sensitive to small 
changes in reaction rate and hence small changes in pressure and gas velocity. Linear 
volatilization rate constants are sufficient to define the steady state recession rate of silica 
formers in combustion environments when in the paralinear regime. Steady state oxide 
thickness, time constant, and recession rate maps allow prediction of SiC and Si.iN 4 
behavior in combustion environments. These maps also provide a framework for 
understanding the limits of the paralinear model. Additional characterization of 
oxidation/volatilization of silica formers is needed to establish the temperature 
dependence of the process. Finally, more studies are needed to understand the chemical 
degradation mechanisms of silica formers in high pressure and very high gas velocity 
combustion environments. 
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Table 1. 

Experimentally determined pressure, temperature, and gas velocity dependence of 
oxidation and volatilization rate constants for SiC in water vapor. 


n, P(H 2 0) n 

T dependence 
(kJ/mol) 

m, v m 

Comments, 

reference 

Parabolic oxidation rate, k 

P 

0.67±0. 19 
0.85±0.05 

28 to 156 

n/a 

not corrected for 0 2 
contribution, (4) 

1.0 

68 

n/a 

silicon, (12) 

0.91 ±0.10 

- 

n/a 

(13) 

— 

35± 1 47 

n/a 

(4) 

Linear vo 

atilization rate, 1 

<1 

1.50±0. 13 

108±7 

0.50±0. 16 

(7) 

— 

56.7±1.7 

— 

(16) 

— 

3+75 

~ 

(4) 


Table 2. 

Summary of relationships describing steady state paralinear oxidation/volatilization. 



Limiting oxide thickness 

Time to steady state 

Steady state 
recession rate 

As a 

function of 
rate 

constants 

x L =k p /2k, 

t L =k p /2kf 

Yl 06 k| 

As a 

function of 
variables 

x L oc(Pv ) _1/2 

t L -(P 2 v) _1 

y L x P 3/2 v 1/2 

Temperature 

dependence 

.. c ..±AH vol -AH ra 'l 

t - rvrl VOl ^ 

v, =exnf“ AHvo1 ) 

L ( RT J 

i t exp 

^ RT ) 

y l ex P 

A R1 L 
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Table 3. 

Summary of test conditions for SiC and Si,N 4 exposed in combustion environments. 


Exposure, 

reference 

material 

P, atm 

P(H 2 0), 

atm 

v, m/s 

T, °C 

Approximate 
time (h) 

A. furnace 1 (3) 

SiC 

1 

0.1 

0.0044 

1200-1400 

100 

B. furnace 2 (4) 

SiC 

1 

0.5 

0.044 

1200-1400 

100 

C. high pressure 
furnace (5) 

SiC 

Si,N, 

10 

1.5 

5x10 4 

1200 

500 

D. mach 0.3 
burner rig (6) 

SiC 

1 

0.1 

100 

1200-1400 

3500 

E. high pressure 
burner rig (7) 

SiC 

5.7 

0.7 

20 

1200-1450 

100 

F. turbine (10) 

Si,N, 

8.7 

0.87 

573 

1066-1260 

815 


Table 4. 

A comparison of measured oxidation and volatilization rate constants 
with rate constants calculated from Equations 10 and 15. 



k„, u 2 SiO,/h 

k„ |U. SiO,/h 

calculated 

measured 

calculated 

measured 

A. laboratory furnace 1 
(3) 

0.044 

0.062 

2.4x1 O' 4 

NA 

B. laboratory furnace 2 
(4) 

0.22 

0.40 

1.9x1 0‘ 2 

1.3x1 0' 2 

C. high pressure 
furnace (5) 

0.66 

NA* 

5.7x1 O' 5 

NA* 

D. mach 0.3 burner rig 
(unpublished work) 

0.044 

NA 

3.6x1 O' 2 

2.0x1 0‘ 2 

E. high pressure burner 
rig (13,7 )t 

0.31 

0.31 

3.3x10' 

3.3x10 1 

F. turbine (10) 

0.39 

NA 

2.2 

0.2- 1.5 


*Under these conditions a non -protective scale is formed on SiC, paralinear 
oxidation/volatilization does not occur and the rate constants are not physically 
meaningful in the context of the paralinear mechanism. 

■(•Calculated and measured values are identical for this case since Equations 10 and 15 are 
based on these measurements. 
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Table 5. 

Summary of calculated and observed oxide thickness for SiC and Si,N 4 exposed in 

combustion environments. 


Exposure, reference 

x L , ft 

t L , h 

x measured, p 

t measured, h 

t>t L 

A. furnace 1 (3) 

91 

4x10 s 

1. 9-2.8 

100 


B. furnace 2 (4) 

6 

300 

2. 5-4.5 

100 


C. high pressure furnace* (5) 

58 

1 x 1 0 4 

25-30 

50-60 

100 

500 


D. mach 0.3 burner rig 
(unpublished work) 

0.6 

17 

2-7 

200 


E. high pressure burner rig (7) 

0.6 

3 

2-6 

100 


F. turbine (10) 

0.1 

2 min. 

none detected 

815 

V 


*Paralinear behavior not observed under these conditions. x L and t L would be meaningful 
if a dense scale was formed. 


Table 6. 

Steady state recession rates calculated from Equation 21 for SiC at 13 16°C under 
conditions in Table 3 assuming paralinear oxidation/volatilization. 



y , p SiC/h 

A. Furnace 1 

2x 1 0 4 

B. Furnace 2 

2x1 0' 2 

C. High pressure furnace 

6x1 0' 3 

D. Mach 0.3 burner rig 

4x1 O' 2 

E. High pressure burner rig 

3x10-' 

F. turbine 

2 


Table 7. 

Limits for paralinear oxidation/volatilization. 


limit 

velocity 

pressure 

controlling recession 
mechanism 

1. negligible 

volatilization/ 

recession 

low 

low 

parabolic oxidation rate 

2. long time 
constant 

low 

low 

parabolic oxidation rate 

3a. non- 
protective scale 

low 

high 

rate of oxidant transport 
through porous scale 

3b. thick 
spalling scale 

low 

high 

parabolic oxidation rate, 
spallation rate? 

4. thin or 
missing scale 

high 

high 

linear oxide growth rate? 
active hydroxide formation? 
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Figure 1. Dimensional change for SiC due to paralinear oxidation and volatilization. 
Oxide thickness and SiC recession curves calculated for k p =0.25 |T/h and k=0.21 |i/h.. 

Figure 2. Determination of the time constant for achieving steady state oxide thickness. 
The solid curve represents paralinear oxide growth. The dotted lines represent oxide 
thickness for parabolic kinetics and the steady state oxide thickness. Oxide thickness and 
time are normalized by x L and t co , respectively. 

Figure 3. Determination of the time constant for achieving the steady state recession rate. 
The solid curve represents paralinear recession. The dotted lines represent recession for 
parabolic kinetics and the steady state recession. Recession depth and time are 
normalized by y n =k p Vk,’ and t cr , respectively. 

Figure 4. Limiting oxide thickness map calculated in terms of the oxidation and 
volatilization rate constants for SiC at 13 16°C. 

Figure 5. Steady state time constant map calculated in terms of the oxidation and 
volatilization rate constants for SiC at 1316°C. 

Figure 6. Limiting oxide thickness map calculated in terms of the total pressure and gas 
velocity for SiC at 1316°C in a hydrocarbon combustion environment. 

Figure 7. Steady state time constant map calculated in terms of the total pressure and gas 
velocity for SiC at 1316°C in a hydrocarbon combustion environment. 

Figure 8. Steady state recession rate map calculated in terms of the total pressure and gas 
velocity for SiC at 1316°C in a hydrocarbon combustion environment. 

Figure 9. Representative oxide cross-sections formed on silica-formers for each of the 
case studies. A) laboratory furnace 1: CVD SiC, 1300°C, 1 atm, 10% H2O/90% O 2 , 0.4 
cm/s, 1 13h. B) laboratory furnace 2: CVD SiC, 1300°C, 1 atm, 50% H2O/50% O 2 , 4.4 
cm/s, lOOh. C) high pressure furnace: CVD SiC (left), AS800 Si 3 N 4 (right), 1200°C, 10 
atm, 15% H 2 O, 0.05 cm/s, 500h. D) mach 0.3 burner rig: sintered SiC, 1300°C, 1 atm, 
10% H 2 O, 100 m/s, 200h. E) high pressure burner rig: CVD SiC, 1316°C, 5.5 atm, 10% 
H 2 0, 20 m/s, 98h. F) turbine: AS800 Si 3 N 4 , 1066-1260°C, 8.9 atm, 10% H 2 0, 162-573 
m/s, 815h. 

Figure 10. Limits for the paralinear oxidation and volatilization model mapped in terms 
of total pressure and gas velocity for SiC at 1316°C in a hydrocarbon combustion 
environment. The dots represent the case studies referred to in Table 3. 
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